Elliptic flow in y/s = 200 GeV Au+Au collisions and y/s = 2.76 TeV Pb+Pb collisions: 
insights from viscous hydrodynamics + hadron cascade hybrid model 
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Using the newly developed hybrid model VISHNU which connects viscous hydrodynamics with a 
hadron cascade model, we study the differential and integrated elliptic flow V2 at different centrality 
bins for 200 A GeV Au+Au collisions and 2.76 A TeV Pb+Pb collisions. We find that the average 
Quark Gluon Plasma (QGP) specific shear viscosity jj/s slightly increases from Relativistic Heavy 
Ion Collider (RHIC) to Large Hadron Collider (LHC) energies. However, a further study assuming 
different temperature dependencies for (^/s)qgp shows that one cannot uniquely constrain the form 
of (v/ s )qgp(T) by fitting the spectra and V2 alone. Based on our current understanding, the question 
whether the QGP fluid is more viscous or more perfect in the temperature regime reached by LHC 
energies is still open. 

PACS numbers: 25.75.-q, 12.38. Mh, 25.75.Ld, 24.10.Nz 
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The first heavy-ion data from the Large Hadron Col- 
lider (LHC) have revealed many phenomena that are very 
similar to those seen at the Relativistic Heavy-Ion Col- 
lider (RHIC) at lower beam energies, such as elliptic flow 
and jet energy-loss. The question arises whether the in- 
creased reach in temperature and energy-density attain- 
able at the LHC allows for the identification of systematic 
trends, for example in the temperature dependence of the 
specific shear viscosity rj/s of the produced Quark-Gluon 
Plasma (QGP). 

Measurements by the ALICE collaboration for 
2.76 A TeV Pb+Pb collisions shows that the charged 
hadron multiplicity density is about a factor of 2.2 higher 
than the one from 200 A GeV Au+Au collisions [J 0. 
Assuming a similar QGP thcrmalization time at RHIC 
and LHC energies and a linear relationship between the 
final multiplicity and initial entropy, one finds that the 
initial temperature of the QGP fireball at lower LHC en- 
ergies is about 30% larger than the one at top RHIC 
energies. Meanwhile, the differential elliptic flow V2(pt) 
of charged hadrons measured by the ALICE collabora- 
tion as a function of transverse momentum px (using 
the 4-particle cumulant method) is, up to transverse mo- 
menta of 3 GeV/c, nearly identical to that measured by 
the STAR collaboration at RHIC, independent of colli- 
sion centrality When integrated over the transverse 
momentum, on the other hand, the total V2 from ALICE 
is about 30% higher than that from STAR, due to an 
increase of the mean px for the LHC spectra Q . 

The larger integrated charged hadron elliptic flow at 
the LHC implies a larger total momentum anisotropy 
H,Q and suggests a higher efficiency of the QGP fluid for 
converting the initial spatial deformation of the collision 
fireball into anisotropic collective flow. This raises the 
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following questions: is the QGP fluid still strongly cou- 
pled at LHC energies? Does the QGP fluid have a similar 
specific shear viscosity at RHIC and LHC energies? In 
this article, we shall address the latter question. 

Several groups have recently published analy- 

ses of the ALICE data Q based on a purely hydrody- 
namic approach, and another group offered an assess- 
ment of the same data within an ideal hydrodynamic + 
hadron cascade hybrid approach [9|. In contrast, we here 
use the newly developed hybrid approach VISHNU [l(| 
which correctly describes viscous and other dissipative 
effects in both the early QGP and late hadronic rescat- 
tcring stages, including the breakdown of chemical equi- 
librium in the hadron gas Q, by connecting the (2+1)- 
dimcnsional viscous hydrodynamic code VISH2+1 
with the microscopic hadronic transport model UrQMD 
[12 (see Refs. f]~0l - |l2| for details). Comparing to a pure 
hydrodynamic approach the kinetic and chemical 

freeze-out of the system is naturally described by the mi- 
croscopic hadronic evolution part of VISHNU through 
elastic, semi-elastic and inelastic scatterings, which elim- 
inates the corresponding freeze-out parameters used in 
purely hydrodynamic calculations. The transition from 
the viscous hydrodynamic to the microscopic evolution 
parts of the model occurs at the switching temperature, 
T sw , which is chosen to be 165 MeV, i.e. near the QGP 
phase transition temperature T c and adjusted to repro- 
duce the chemical freeze-out temperature at RHIC ener- 
gies [l3|. This is the highest T sw possible for a hadronic 
transport description as well as the lowest T sw possible for 
hydrodynamics without introducing additional parame- 
ters related to hadronic non-chemical equilibrium and 
dissipative effects in the hydrodynamic part [Io| . For the 
hydrodynamic evolution above T sw , we use the equation 
of state (EOS) s95p-PCE Q which is based on recent 
lattice results [lj|. 

In this manuscript, we will focus on the specific shear 
viscosity of the QGP at RHIC and LHC energies, ne- 
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FIG. 1: (Color online) Left: pr spectra for all charged hadrons. Experimental data are from ALICE [33 
Pt spectra for pions and protons. Experimental data are from STAR [34[ and PHENIX [35j . 



and STAR [32]. Right: 



glecting the bulk viscosity 1 and assuming zero net baryon 
density and heat conductivity. In our calculations we 
shall either use a constant r//s at RHIC and LHC ener- 
gies (which represents the averaged shear viscosity effects 
during the whole evolution of the QGP fireball) or input 
different temperature dependencies for (r\j s)qgp(T) in 
the temperature range covered by both RHIC and LHC 
energies. The corresponding relaxation time is set as 
TV = 3 V /(sT) 2 

In order to start a VISHNU calculation, initial condi- 
tions are required. The two most popular geometric 
models for initial particle production including fluctu- 
ation effects for high-energy heavy ion collisions are the 
Monte-Carlo Glauber Model (MC-Glauber)JpJl|| and 
the Monte-Carlo KLN Model (MC-KLN) {Ulllf Fol- 
lowing Ref . [1, H, , we account for event-by-event fluc- 
tuations on average by using an initial entropy density 
profile that is the result of averaging over a large num- 
ber of fluctuating initial entropy density distributions. 
Such an average can be done in two ways: either by re- 
centering and rotating each Monte-Carlo event to align 
the major and minor axes of each initial density distri- 
bution before averaging (initialization in the participant 
plane) or by directly averaging without re-centering and 
rotating (initialization in the reaction plane). In this pa- 
per, we use the reaction plane method since we will com- 
pare our theoretical results with the elliptic flow data 
extracted by STAR and ALICE using the 4-particle cu- 
mulant method, «2{4}, which minimizes non-flow effects 
[2l| and measures t>2 in the reaction plane under the as- 
sumption of Gaussian fluctuations |22j . The assumption 



of Gaussian fluctuations was challenged in Ref. [23( for 
the most central and peripheral collisions; a recent analy- 
sis [24| confirmed, however, that for collision centralities 
up to about 40% reaction-plane averaged initial condi- 
tions produce initial eccentricities that are very close to 
e{4}, the presumed driver for i>2{4} [22| . Only for more 
peripheral collisions does this shortcut break down (2~4| . 

Questions about the validity of the direct compar- 
ison of the reaction-plane averaged theoretical results 
with ^{4} data from peripheral collisions engender cor- 
responding uncertainties for the extracted values of the 
QGP viscosity. Additional 5—10% uncertainties arise 
from the recent observation (23 - [26j (see also (27}) that 
single-shot hydrodynamic evolution of a smooth aver- 
aged initial profile (as employed here) slightly overpre- 
dicts the elliptic flow compared to an event- by- event hy- 
drodynamic evolution of each fluctuating and highly in- 
homogeneous initial profile separately where the average 
over the event ensemble is taken only at the end. In con- 
trast to Ref. [1, [H, we here do not aim to extract the 
QGP viscosity at RHIC and LHC energies with reliable 
uncertainty estimates, but rather to investigate the rel- 
ative change of the QGP viscosity from RHIC to LHC 
energies. Assuming that distortions from non-Gaussian 
eccentricity fluctuations and the replacement of event- by- 
event evolution by single-shot hydrodynamics are similar 
for the STAR and ALICE data, their presence will not 
affect the conclusions regarding the variation of the QGP 
viscosity with collision energy drawn from a comparison 
of i>2{4} data with our calculations. For the same reason 
and for the sake of computing efficiency, we only utilize 
the MC-KLN initialization in the reaction plane, 3 and 



1 Ref. flEll shows that the bulk viscosity results in less than 20% 
contamination for the extracted QGP shear viscosity value, due 
to the critical slowing down of the bulk relaxation time near T c . 

2 The choice of tv in the QGP phase has small to negligible effects 
on the final spectra and elliptic flow [TUl . 



3 Studies show that the choice of MC-Glauber or MC-KLN ini- 
tializations will only affect the absolute value of the extracted 
specific shear viscosity at RHIC and LHC energies, but that the 
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leave the case for a detailed extraction of QGP viscosity 
at LHC energies to future studies. 
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FIG. 2: (Color online) Centrality dependence of the 
charged hadron pseudo-rapidity density per participant pair 
(dN/drf)/(N Dart /2). Experimental data are from ALICE 0|, 
STAR [3(| and PHOBOS [H. Theoretical lines in both pan- 
els are from VISHNU with different constant rj/s as input (see 
the text for the details of other inputs and parameters). 



The initial time To and the normalization of the aver- 
aged initial entropy density profile need to be fixed from 
experimental data. Following Refs. [1, H|, we use the fol- 
lowing parameter sets for the shear viscosity to entropy 
ratio rj/s and hydrodynamic starting time to: (0.16, 0.9 
fm/c), (0.20, 1.05 fm/c) and (0.24, 1.2 fm/c). Please note 
that for a larger value of the QGP viscosity, we use a later 
starting time to to compensate for the additional radial 
flow generated by that larger viscosity [1, Hj| . After tun- 
ing the normalization of the initial entropy density to 
approximately reproduce the final state charged hadron 
multiplicity per unit of pseudo-rapidity in 200 A GeV 
central Au+Au collisions {dN/drj ~ 690 USE!) and in 
2.76 A TeV central Pb+Pb collisions (dN/drj ~ 1600 Q), 
we find that our calculation provides a good description 
of the data on px spectra for all charged hadrons in most 
central collisions for STAR [H HI and ALICE [H as 
shown in Fig. QJi. Ref. Q also shows that with the above 
parameters one can obtain a good fit to the p^-spectra for 
identified hadrons (such as pions and protons) from most 
central collision to most peripheral collisions at RHIC en- 
ergies. We find that, with the above adjustment of the 
starting time To when changing rj/s, these p^-spectra are 
rather insensitive to the QGP viscosity. In Fig. [T|d, we 
show the pT-spectra for pions and protons in most central 
collisions, and compare the RHIC results with the STAR 
[3! and PHENIX data (H. Due to the current lack 



of ALICE data for identified hadrons, the corresponding 
LHC results are predictions. 

In the MC-KLN initialization, we use the standard 
parametrization for the saturation scale Q 2 A as shown 
in (l8l | , which is tuned to reproduce the centrality depen- 
dence of the charged hadron multiplicity for 200 A GeV 
Au+Au collisions. Fig. [2] shows that such aparametriza- 
tion also leads to a good description for the slope of the 
(dN ch /drj)/(N part /2)-N part curve in 2.76 A TeV Pb+Pb 
collisions. However, we have to point out that the same 
parametrization of Q 2 A will lead to a slight overpredic- 
tion of the value of dN^/drj at LHC energies as shown in 
Ref. Q • To avoid the over-generation of elliptic flow from 
over-predicted final multiplicities, we tune the normaliza- 
tion of initial entropy density as described above to fit the 
dN/drj in the 0—5% centrality bin. This leads to a good 
fit on the overall magnitude of the (dN^/ 'drf) / '(A r part /2) 
vs. .ZVpart curve. 

Having fixed all parameters, we calculate the differen- 
tial elliptic flow at RHIC and LHC energies for different 
constant values of rj/s as input. Fig. [3] shows that with 
rj/s = 0.16, VISHNU nicely fits the STAR w 2 (pt){4} data 
from to 2 GeV for different centrality bins. In contrast, 
the same rj/s = 0.16 significantly over-shoots the AL- 
ICE U2{4} data (for pt > 0.5GeV/c). After increasing 
rj/s to 0.20 - 0.24, VISHNU can roughly fit the ALICE 
data at higher pt, but still undcr-prcdicts the data for 
Pt < 0.5GeV/c. This effect of undcr-prediction of the 
low pt data is also found in other hydrodynamics-based 
calculations, including the (3+l)-d ideal hydrodynamics 
+ hadron cascade simulations by Hirano et. al 0, the 
(2+l)-d viscous hydrodynamic calculations by Bozek Q 
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trend on how the QGP viscosity changes from RHIC to LHC is 
very similar for MC-Glauber and MC-KLN initializations [28ll . 
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FIG. 3: (Color online) v 2 (p T ) at 10-20%, 20-30%, 30-40% and 
40-50% centrality. Experimental data are from STAR [2(| 
and ALICE 0] obtained from 4 particle cumulant method. 
Theoretical lines are from VISHNU calculations with different 
constant rj/s as input. See text for details. 
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and the event-by-event simulations with a (3+l)-d vis- 
cous hydrodynamics with fluctuating initial conditions 
done by Schenke et al. 0]. In Rcf. Q, the deviation is 
interpreted to be due to non-thermalized particles stem- 
ming from jet fragmentation. While the origin of this 
deviation is still under debate, we conclude from Fig. [3] 
that one needs a larger averaged QGP specific viscosity 
to fit the ALICE ^(pt) at pr > 0.5 GeV/c than the one 
used to fit the corresponding STAR data. This conclu- 
sion rests on the assumption that the theoretical model 
correctly describes the slopes of the charged hadron px- 
spectra at all the centralities shown in Fig. [3] 

Using VISHNU we find that, even at LHC energies where 
almost all of the final momentum anisotropy is generated 
hydrodynamically in the QGP stage, the charged hadron 
V2(pt) continues to grow somewhat during the hadronic 
stage. This hadronic increase of V2{pt) is smaller at 
the LHC than at RHIC, in agreement with earlier find- 
ings using an ideal hydrodynamic + cascade hybrid 
code. At RHIC energies, some of this hadronic increase 
is driven by the creation of additional overall momen- 
tum anisotropy which has not yet quite saturated in the 
QGP phase. At LHC energies, it is mostly caused by 
a hadronic redistribution of the momentum anisotropy 
already established in the QGP phase in px and among 
the various different hadronic species, due the hadronic 
increase in radial flow that pushes t>2 to larger pt, espe- 
cially for heavy particles [33]. This effect is sensitive to 
the chemical composition in the hadron gas (38U4H , and a 
corresponding hadronic increase of V2{pt) is not observed 
in purely hydrodynamic calculations with an equation 
of state that (incorrectly) assumes chemical equilibrium 
among the hadrons even below the chemical decoupling 
temperature T chom w 165 MeV ||. 

Fig. 2] shows the comparison of experimental and theo- 
retical integrated 1)2, obtained from integrating «2(pt) 
with the corresponding pt spectra as weighting func- 
tions. Following the STAR \2j^ and ALICE analy- 
sis, we use the same px and pseudo-rapidity cut in our 
VISHNU calculations (0.15 < p T < 2GeV/c and \rj\ < 1 
at RHIC energy, and 0.2 < p T < 5 GeV/c and \rj\ < 0.8 
at LHC energy). One finds that VISHNU is capable of 
fitting the experimental data with 77/s = 0.16 at RHIC 
and 77/s = 0.20 at LHC, except for the most peripheral 
centrality bins. Comparing our calculation to the STAR 
V2{pt) data with 77/s = 0.16, the solid purple curve with 
square symbols is slightly above the STAR data due to 
the slight overprediction of the pr-spectra around 1 GeV 
as shown in Fig. QJi. Similarly, the value of 77/s = 0.20 
from the fit to the ALICE integrated V2 is slightly below 
the extracted value of 77/s = 0.22 from ALICE V2(pr), 
mainly because of the under- fitting of the ALICE V2(pr) 
at lower px ■ 

In Ref. Q, we discussed that integrated v% is better 
suited than differential V2 for the extraction of the QGP 
viscosity, due to it being directly related to the fluid mo- 
mentum anisotropy and insensitive to other details of hy- 
drodynamic calculation such as chemical components of 
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FIG. 4: (Color online) Integrated V2 as a function of centrality. 
Same illustrations for theoretical and experiential lines as in 
Fig. El 



the hadronic phase, the form of non-equilibrium distribu- 
tion function 5f, bulk viscosity and so on. However, due 
to the current deviation between calculations and the AL- 
ICE V2(pt) data at lower pr, which translates into cor- 
responding errors for the integrated i>2, any extraction of 
the QGP viscosity from the integrated V2 measurements 
alone at LHC energies cannot be considered robust. 

Although both Figs.[3]and|4]indicate that the averaged 
QGP specific viscosity (constant 77/s) slightly increases 
with collision energies, it has to be pointed out that us- 
ing one constant value of 77/s to fit RHIC data and a 
different constant value of 77/s to fit LHC data is not log- 
ically consistent. In other words, one can not describe the 
QGP fluid created at at RHIC energies with 77/s = 0.16 
(T c < T < 2T C ) and then use 77/s = 0.22 (T C <T< 3T C ) for 
the one created at LHC energies. It is a temperature de- 
pendent rj/s(T) that reflects the intrinsic properties of 
the QGP fluid, and this temperature dependence should 
be unique and describe the data both at RHIC and LHC 
energies. 

However, Fig. [5] shows that one can at least find two 
different functional forms of (77/ s)qqp(T), with which 
VISHNU can simultaneously fit the STAR and ALICE 
v 2 (pt) at 30-40% centrality bins. 4 (r]/s)(T)^ a '> monotoni- 
cally increases with T in the QGP phase while (77/s) (T) ^ 
first increases with T and then decrease with T at even 
higher temperature. 5 Please note that the minimum val- 



4 (t?/.s)(T)M and (»?/s)(T)W can also nicely fit the p T -spectra 
for identified hadrons at 30-40% centrality which, due to lack 
of data, we obtained theoretically from VISHNU, using constant 
77/s as input. With px spectra and i>2 (pt) fitted, one can also 
roughly fit the integrated V2, since the latter is calculated from 
the former two. 

5 For the purpose of demonstration, we have chosen simple (even 
unrealistic) forms for the temperature dependence of {ri/s)(T)( a ^ 
and (r)/s)(T)^ with 2-4 free parameters that can easily be fitted 
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ues of (r]/s)(T) in case (a) and (b) are below 0.16, and 
the the maximum values are well above 0.24. 6 Although 
(r]/s)(T)^ shows a smaller and subsequently negative 
slope at higher temperature, it has a higher minimum 
value of 77/s than {ri/s){T)^ a \ 

The fireball evolution may differ between (77/ s)(T)( a ) 
and (r]/s)(T)( b \ However, both the final and particle 
spectra are sensitive only to the time- integral of the QGP 
evolution, which in both cases apparently is very similar 
to the evolutions at fixed values of 77/s = 0.16 at RHIC 
energy and 77/s = 0.22 at LHC energy. Based on our 
current analysis, utilizing only V2 and particle spectra, 
one cannot unambiguously determine the functional form 
of (ji/s)(T) and whether the QGP fluid is more viscous 
or more perfect at LHC energy. 

Our parametrization (77/ s)(T)^ is similar to one of 
the forms used in whose authors studied the sensi- 
tivity of elliptic flow measurements at RHIC and LHC 
to a temperature-dependent increase of (?7/s)qgp above 
the quark-hadron phase transition. The authors of [42[ 
concluded that elliptic flow measurements at RHIC en- 
ergies are insensitive to such an increase, and provided 
convincing evidence that this conclusion does not de- 
pend on their treatment of the late hadronic stage which 
they evolved hydrodynamically, with sudden Cooper- 
Frye freeze-out, rather than microscopically as we do 
here. In contrast to their work, we do find such a sensi- 
tivity; this is why in (?7/s)(T)( a ) we had to lower the min- 
imal 77/s value at low temperatures below the constant 
value of 0.16 that we used in Fig. [3J in order to pre- 



0.25 




serve the theoretical description of the V2(j>t) measured 
by STAR. We have traced this difference to a sensitivity 
to initial conditions: whereas Niemi et at. [42j initialize 
the shear viscous pressure tensor ■n iiv at either zero or 
a non-zero initial value that is independent of 77/s, we 
use Navier- Stokes initial conditions ir^ = 2 (-| (T)j sa^ 
(a^ is the velocity shear tensor [ll|) where it' 11 ' increases 
with T (i.e. towards the fireball center) if 77/s does so. 
This increases the initial transverse pressure gradients 
and thereby affects the final radial and elliptic flow. 

In summary, we have studied the differential and in- 
tegrated V2 at different centrality bins for 200 A GeV 
Au+Au collisions and 2.76 A TeV Pb+Pb collisions, us- 
ing the hybrid model VISHNU which describes the expan- 
sion of a QGP using viscous hydrodynamics and the suc- 
cessive evolution of hadronic matter with a microscopic 
transport model. We find that, in order to describe the 
STAR and ALICE «2{4} data with reaction-plane av- 
eraged MC-KLN initial conditions, one needs an aver- 
aged QGP viscosity 77/s w 0.16 at RHIC energies and 
77/s ?s 0.20—0.24 at LHC energies. Although this result is 
in qualitative agreement with expectations from weakly 
coupled QGP calculations [43[ and from recent lattice 
simulations |44j, both of which shows that the specific 



QGP shear viscosity increases with temperature, a more 
detailed analysis with VISHNU utilizing a temperature de- 
pendent (t7/s)qgp(T 1 ) as input shows that one cannot 
uniquely constrain the form of (rj/s)(T) by fitting the 
spectra and V2 alone. Based on our phenomenological 
approach, it remains an open question whether the QGP 
fluid is more viscous or more perfect in the temperature 
regime reached by LHC energies. 
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FIG. 5: (Color online) Two example of temperature depen- 
dent (7?/s)qgp(T), with which VISHNU can simultaneously fit 
the STAR and ALICE v 2 (p T ) at 30-40% centrality. 



from spectra and V2 at RHIC and LHC energies. 
Since below T c ^ Bm = 165 MeV the fluid is described microscopi- 
cally, the behaviour of (rj / s)(T)( a ^ shown in Fig. [5] in the re- 
gion T < 165 MeV is irrelevant for our calculation. We made 
no attempt to model in detail the (presently unknown) form of 
{w/s)(T) in the phase transition region. 
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